Theory of Activated Transport in Bilayer Quantum Hall Systems 
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We analyze the transport properties of bilayer quantum Hall systems at total filling factor u — 1 in 
drag geometries as a function of interlayer bias, in the limit where the disorder is sufficiently strong 
to unbind meron-antimeron pairs, the charged topological defects of the system. We compute the 
typical energy barrier for these objects to cross incompressible regions within the disordered system 
using a Hartree-Fock approach, and show how this leads to multiple activation energies when the 
system is biased. We then demonstrate using a bosonic Chern-Simons theory that in drag geometries, 
current in a single layer directly leads to forces on only two of the four types of merons, inducing 
dissipation only in the drive layer. Dissipation in the drag layer results from interactions among 
the merons, resulting in very different temperature dependences for the drag and drive layers, in 
qualitative agreement with experiment. 

PACS numbers: 73.43-f, 03.75.Lm, 73.21-b 



Double layer quantum Hall systems at filling factor 
V = 1 display many properties akin to those of superflu- 
ids . This behavior results from the pairing of electrons 
in one layer with holes in the other, producing excitons 
that condense into a state with interlayer coherence even 
in the absence of tunneling . In experiments these sys- 
tems display a strong interlayer tunneling peak at zero 
bias, reminiscent of the DC-Josephson effect Q, and van- 
ishing single-layer resistances as temperature T — > in 
counterflow experiments Nevertheless, the "superflu- 
idity" in this system remains imperfect: there is no truly 
dissipationless transport at low but finite temperature in 
either of these types of experiments. Recently, it has been 
demonstrated that this behavior may be qualitatively un- 
derstood if one assumes that disorder produces unpaired 
merons - the analog of vortices in a thin film superfluid - 
that remain weakly mobile at any finite temperature Q . 

While these experiments strongly suggest the near- 
coherence of the two layers in this system, one class of 
experiments has so far defied explanation in these terms. 
These are drag measurements, in which current is in- 
jected and removed in a single layer, and the voltage drop 
measured in either layer. The resulting resistances when 
measured as a function of temperature have roughly ac- 
tivated behaviors. The activation energies for the drive 
and drag layers behave very differently with interlayer 
bias: the former are asymmetric with respect to the bias 
direction, while the latter are roughly symmetric 0, 01 ■ 
Naively this suggests that each layer has separate quasi- 
particles with different activation energies, and interlayer 
coherence essentially plays no role. Yet such a picture is 
very difficult to reconcile with the experiments described 
above, in which imperfect superfluidity is manifest. 

In this paper, we will propose a solution to this puzzle. 
Our approach involves a transport theory for this system 



in which disorder is incorporated Q via a slowly-varying 
random potential (such as results from a remote dop- 
ing layer), producing puddles of charged quasiparticles 
Q. In the context of quantum Hall bilayers, one expects 
such quasiparticles to be constructed from topological 
defects These merons carry a vorticity (s = ±1) 
in the relative phase of the wavefunction amplitude for 
each layer, and an interlayer polarization (a = ±1) aris- 
ing from a charge imbalance in the meron core that may 
tilt into either of the two layers. Moreover, because of 
the remarkable properties of electrons in a single Landau 
level [9], the topological charge of such objects is tied to 
their physical charge: q^a =^ ^sav^^, with v^i — vu 
and v-i = ul, and i^u(L) the filling factor in the upper 
(lower) layer. This will have important consequences for 
their response to currents. 

In this "coherence network" model most of the 
area is taken up by puddles where both the interlayer 
coherence and the incompressibility are lost due to the 
large local density of merons. However, puddles are sep- 
arated by narrow strips of incompressible Hall fluid [8| 
with local filling factor very close to 1. Thus, any mea- 
sured activation energy actually refiects the energy bar- 
rier for a meron to cross such an incompressible strip 
in moving from one puddle to its neighbor. Below, we 
describe Hartree-Fock (HF) calculations to estimate the 
dependence of this energy barrier on interlayer bias. Our 
results demonstrate that the activation barrier depends 
significantly on the relative orientation of the polariza- 
tion of merons and the bias. This explains the multiple 
activation energies in quantum Hall bilayers. 

To understand the transport properties of the merons, 
we map the system to a bosonic superfluid state using a 
Chern-Simons (CS) transformation, in which each elec- 
tron is understood as a boson with a single unit of mag- 
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netic flux directed opposite to the applied magnetic field 
At mean- field level the magnetic field is can- 
celed, and the quantum Hall state may be understood 
as a Bose condensate of the composite bosons. In such 
a description, the quantum Hall bilayer is condensed in 
two senses: with respect to the (bosonic) charge degrees 
of freedom, and with respect to the interlayer degree of 
freedom. In computing the total force on a meron due 
to currents in the system, one must account for forces 
due to the excess charge on a meron, with which there 
is an associated magnetic flux, and due to its vorticity 
in the interlayer phase. In the case of current in just 
one layer these forces cancel precisely for merons with 
polarization directed toward the current-carrying layer, 
so that only half the merons move in direct response to 
the current. Moreover, the induced voltage due to mo- 
tion of the driven merons turns out to lie completely in 
the driven layer. Then the induced voltage depends on 
activation barriers for merons of a single polarization, 
leading to an asymmetric activation energy with respect 
to bias, as seen in experiment. 

A voltage drop is induced in the drag layer only 
through interactions: a driven meron may pair with 
an undriven meron of opposite polarization (to form a 
bimeron [l|) over some distance, inducing a voltage in 
the drag layer. Since the relevant activation energy now 
involves the crossing of merons of both polarizations over 
incompressible strips, one expects the resulting activa- 
tion energy to be symmetric with bias. This again is the 
experimental observation. 

Numerical calculation of energies: In equilibrium, one 
expects the energy of merons residing in different puddles 
to be roughly equal. When a meron hops from one puddle 
to the next it must pass through an incompressible strip, 
where its energy is higher. Computing the activation 
energy directly for such a process is challenging because 
it is difhcult to fully model the effects of the puddles. In 
what follows we will estimate the electrostatic (Hartree) 
contribution to the activation energy, and demonstrate 
that it is sensitive to the interlayer bias. 

Our approach is a Green's function equation of motion 



method 11| whose application to meron-antimeron 



states of a quantum Hall bilayer has been described else- 
where UU- The method generates order parameters 
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where c -^ ^ creates an electron in layer i in a lowest 
Landau level (LLL) state with guiding center quantum 
number X, for Hartree- Fock states with crystalline 
order, characterized by a set of reciprocal lattice vectors 
{G}. To inject merons and antimerons into this state, 
one works slightly above or below filling factor 1. (We 
use V = 1.02 in the results described here). In the limit 
of very small tunneling, the lowest energy state of the 
system is a square lattice, with merons at the center 
and corner of the unit cells, and antimerons at the face 




FIG. 1: The z-component of pseudospin density (in units of 
l/27r^^) in a square lattice of merons with two electrons per 
unit cell at d/i = 1.0 and u — 1.02. Black lines are schematic 
depiction of the imposed barrier potentials. 



centers. The results can be conveniently expressed in 
terms of a pseudospin vector S, defined by S^+iSy = pi2 
and = pn — P22- In this language a meron lattice is 
a complicated non-coUinear ferromagnetic state. Fig. [1] 
illustrates 5*2 for a typical such lattice. Note that the 
polarizations of the merons, represented by the sign of 
Sz at their cores, is opposite for merons and antimerons, 
so that the charge of the two objects is the same. 

While the puddles of merons largely screen the local 
disorder potential, inside the incompressible strips there 
is no such screening. Moreover, because of their high den- 
sity of merons, the pseudospin stiffness inside the pud- 
dles is significantly compromised relative to that of the 
incompressible strips separating them. Thus the maxi- 
mum energy configuration for a meron crossing a strip 
will occur when the meron is centered upon it. (Note 
that the strip widths are of order ^ [1, [1] , which is nar- 
row compared to the meron size.) Our HF approach al- 
lows us to investigate the effect of interlayer bias on the 
first of these contributions to the energy barrier. To do 
this, we add an external "box" potential of height Vq 
and width £ along finite strips in the unit cell, forming a 
checkerboard pattern (see Fig. [T]). These potential strips 
represent the difference in electrostatic energy for charge 
located in a strip and charge located in a puddle. By 
adjusting the position of the meron lattice relative to the 
external potential, we can arrange for the potential to be 
directly beneath the a = +1 merons, the a = —I merons, 
or between merons. The energy difference between the 
"on-meron" and "off-meron" configurations yields an es- 
timate of the electrostatic energy barrier for a meron of 
a particular polarization a to cross a strip. In practice, 
we adjust the barrier height Vo so that this energy differ- 
ence matches the activation energy in Ref. ^ for unbiased 
wells. We can then examine how this activation energy 
evolves as the layers are biased. 

Fig. [5] illustrates some typical results. Merons with 
polarization oriented in the direction opposite that of 
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the bias have an increasing activation energy with bias, 
while those with polarization in the same direction de- 
crease in energy. (Data for Av < were generated us- 
ing Av > results, which is valid due to a symmetry 
upon interchange of the layers.) The result may be un- 
derstood qualitatively by noticing that the charge of the 
former (latter) increases (decreases) with bias, leading 
to a larger (smaller) electrostatic energy cost for travers- 
ing the barrier. Thus, as observed in experiment and has 
been so challenging to explain, one obtains two activation 
energies with opposite slopes as a function of Av. 

The left inset of Fig. [2] illustrates the slope of such 
curves at Av = for different layer separations. The 
apparent non-monotonic behavior may be understood as 
follows. As d increases from small values, the merons 
become smaller [ij so that a larger portion of their area 
lies in the barrier region as they cross. When d becomes 
still larger, the interlayer exchange interaction decreases 
enough that the meron charge density deforms when cen- 
tered on the barrier. With less charge in the barrier re- 
gion, the sensitivity to bias decreases. This this behavior 
is a testable prediction of our model. 

Finally, the right inset of Fig. [5] illustrates the activa- 
tion energy at zero bias for fixed d and various values of 
field B {£ (X yJl/B). One may see that the activation 
energy in our model decreases with increasing B, as is 
the general trend in experiment. 

The slope of the activation energy as a function of bias 
in our model is smaller than that found in experiment Q , 
leading to an activation energy about a factor of 2 smaller 
than what is observed experimentally at 10% polariza- 
tion. Although there may be several reasons for such 
a discrepancy, a large component of it is likely to come 
from the absence of quantum fluctuations in our model. 
This compromises the pseudospin stiffness in the incom- 
pressible region, and since the experiments are operated 
near the coherent-incoherent transition for the system 
0, the effect should be considerable. Application of bias 
is known to increase coherence effects in bilayers 
implying that the stiffness in the barrier region should 
significantly increase with bias. This will increase the 
sensitivity of the activation energy barrier to the bias. 

Meron Dynamics: Having seen how merons may dis- 
play multiple activation energies when the bilayer is bi- 
ased, we now turn to a description of their dynamics and 
the associated dissipation. We work in the composite bo- 
son picture, in which the electron system at v = l\s mod- 
eled as a Bose condensate in zero average magnetic field, 
with an additional degree of freedom, the pseudospin. 

In general, any uniform current in a bilayer system 
can be decomposed into counterflow (CF) and co-flow 
currents. Using the superfluid analogy, CF may be de- 
scribed by spatial rotation of the order parameter phase, 
which in the pseudospin description is the argument of 
Sx + iSy. This same phase angle may alternatively be 
interpreted as the condensate phase when the interlayer 



coherent state is described as an exciton condensate |2f. 
Because of the vorticity associated with a meron, such 
currents induce Magnus forces. In addition, since merons 
are charged objects, in the composite boson picture they 
carry residual magnetic flux. Coflow current therefore 
induces a Lorentz force on a meron in direct analogy to 
flux lines in superconductors [isj . A non-zero net force 
on a meron will induce motion, which in turn creates a 
voltage in each layer. Since the effective excess flux car- 
ried by a meron is q^o, with $o = he/ e, the Lorentz force 
from coflow current may be written [15l | as 

F = {Uo)Jxz, (1) 
c 

in which J is the sum of upper and lower layer current 
density, Ju + Jl- To compute the CF force, we need to 
know the effective velocity of the condensed particle-hole 
pairs, relative to the average velocity of all the electrons 
[lO|, Vs = 2ttP{Ju + Jl). Without loss of generality, one 
may assume < vjj so that the exciton velocity in the 
lab frame, Vexi is the electron velocity in the lower layer. 
Then in the frame comoving with the average electron 
velocity, the CF current density is J(?p = ^^{vex — 
Vs) = Jl — i^l{Jl + Ju)- Since the meron contains a full 
unit of vorticity of the excitonic condensate phase, the 
force due to CF current becomes 

FcF = hsJ^°P xz = {es'pQ/c)J^°P x z (2) 

in which s is the vorticity. Adding ([1]) and ^ and using 
q = —sav^cr yields the total force on a meron [16|, 

-* es -* — ' 

^T = y$o[(l + a)JL-(l-(T)Jc/] xz. (3) 

From this expression, it is immediately apparent that 
merons of a given polarization cr = ±1 respond only to 
the current in a single layer. 

Connection to Experiment: The force Fs^o on merons 
of vorticity s and polarization a will cause them to flow 
with a velocity Ug^a = l^-s.aFs^a where fj,s.ij is an effective 
mobility, which we expect to be thermally activated, with 
a bias dependence of the activation energy as discussed 
above. We use the resulting motion of the vortices to 
find the voltage drops between two points a distance yo 
apart along the direction of electron current. From the 
Josephson relation, due to the vorticity in the interlayer 
phase an interlayer voltage drop AV is induced when 
merons pass between these two points [5], 

AV^AVu-AVL = -—yoy"n,^su,^, (4) 

where Us^^r is the meron density. On the other hand, the 
CS magnetic flux moving with the merons between the 
two points induces a potential drop between electrons 
near the two points that is independent of the layer in 
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which they reside. This leads to the condition 



--yo > 



(5) 



In a drag geometry we have, for example, Jl = and 
Jjj = with / the total current and W the sample 
width. Combining Eqs. [HandO we obtain AVl = and 



AVu 



yo 



(6) 



/ W 

Notice the final result depends on the mobility of only 
merons with polarization cr = — 1. It immediately follows 
that the voltage drop in the drive layer is asymmetric 
with respect to bias, precisely as observed in experiment. 

In order to explain the voltage drop in the drag layer 
{AVl 7^ 0) we must identify how forces on the a = +1 
merons might arise. A natural candidate for this is the 
attractive interaction between merons with opposite vor- 
ticities, which in the absence of disorder binds them 
into pairs at low meron densities. Assuming that driven 
merons crossing incompressible strips will occasionally 
be a component of these bimerons, a voltage drop in the 
drag layer will result. The mobility of such bimerons is 
limited by the energy barrier to cross an incompressible 
strip. Since these strips are narrow compared to the size 
scale of the constituents of the bimeron, we expect the 
activation energy to be given approximately by the max- 
imum of the activation energies for merons of the two 
polarizations cr = ±1. This leads to a drag resistance 
much smaller than that of the drive layer, and an activa- 
tion energy that is symmetric with respect to bias. Both 
these behaviors are observed in experiment [l8| . 

It is interesting that, within our model, as sam- 
ples become increasingly clean one expects such mobile 
bimerons to become more prominent relative to single 
merons, so that the voltage drop in the two layers will in- 
creasingly match with decreasing disorder and decreasing 
temperature. In principle a drag experiment in a sample 
clean and cold enough that all merons are paired will re- 
sult in precisely the same voltage drop in each layer, so 
that pure countcrflow becomes dissipationless, and coflow 
dissipation occurs in a manner such that one cannot dis- 
tinguish whether the electrons are moving in the upper 
or lower layer. Similar behavior should ensue with in- 
creased interlayer tunneling, for which merons and an- 
timerons become more strongly bound into pairs, forming 
the quasiparticles of the system (bimerons) [H, 0, 13]. 
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FIG. 2: Activation barrier for merons as a function of layer 
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vs. d at zero bias. Right inset: Activation energy vs. d/£ at 
zero bias, for fixed d {d = t sX B = IT). 
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